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Abstract. A review of meson emission in heavy ion collisions at incident 
energies from SIS up to collider energies is presented. A statistical model 
assuming chemical equilibrium and local strangeness conservation (i.e. strangeness 
conservation per collision) explains most of the observed features. 

Emphasis is put onto the study of K + and K~ emission at low incident 
energies. In the framework of this statistical model it is shown that the 
experimentally observed equality of K + and K~ rates at "threshold-corrected" 
energies \fs — ^/sth ls due to a crossing of two excitation functions. Furthermore, 
the independence of the K + to K~ ratio on the number of participating nucleons 
observed between SIS and RHIC is consistent with this model. 

It is demonstrated that the K~ production at SIS energies occurs 
predominantly via strangeness exchange and this channel is approaching chemical 
equilibrium. The observed maximum in the K + /7r+ excitation function is also 
seen in the ratio of strange to non-strange particle production. The appearance of 
this maximum around 30 A-GeV is due to the energy dependence of the chemical 
freeze-out parameters T and n b ■ 



1. Introduction 

Central heavy ion collisions at relativistic incident energies represent an ideal tool to 
study nuclear matter at high temperatures. Particle production is - at all incident 
energies - a key quantity to extract information on the properties of nuclear matter 
under these extreme conditions. Particles carrying strangeness have turned out to be 
very valuable messengers. 

A specific purpose of this paper is the presentation of the evolution of strange 
particle production over a large range of incident energies. The data at low incident 
energies, i.e. at and below the production threshold in NN collisions will be shown 
in more detail as these results are less presented in other works. Many results are 
shown together with a theoretical interpretation. The attempts to describe particle 
production yields with statistical models |2[ [| ||, [|, |[ 0, |§| have turned out to be 
very successful over this large domain of incident energies. 
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2. General Trends 



2.1. Production of K+ and K~ from SIS to RHIC 

At incident energies around 1 A-GeSf pion and kaon production is very different: 
Pions can be produced by direct NN collisions but kaons not. The threshold for K + 
production in NN collisions is 1.58 A-GeV and only collective effects can accumulate 
the energy needed to produce a K + together with a A (or another strange particle) 
for strangeness conservation. 

The measured multiplicities for pions and K + differ evidently by orders of 
magnitude. They exhibit a further very pronounced contrast: While the pion 
multiplicity per number of participating nucleons A„ art remains constant with A part , 
the K + multiplicity per A part rises strongly (Fig. [y). The latter observation seems 
to be in conflict with a thermal interpretation, which - in a naive view - should give 
multiplicities per mass number A being constant. 

Usually, the particle number densities or the multiplicities per A part , here for 
pions, are described in a simplified way by a Boltzmann factor 

M n ( <E 7I > 

exp 



Apart \ T 

with the temperature T and the total energy < >. 

The production of strange particles has to fulfil strangeness conservation. The 
attempt to describe the measured particle ratios including strange hadrons at AGS and 
SPS using a strangeness chemical potential fj,$ is quite successful (I], ||, ^| p|. However, 
this grand-canonical treatement is not sufficient, if the number of produced strange 
particles is small. Then, a statistical model has to take care of exact strangeness 
conservation in each reaction as introduced in fic| . This is done by taking into 
account that e.g. together with each K + a A or another strange particle is produced: 
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where T is the temperature, fj,B the baryo-chemical potential, gi the degeneracy 
factors, V the production volume for making the associate pair (see j|, ||) and Ei 
the total energies. We note that this volume is not identical to the volume of the 
system at freeze out. The volume parameter V is taken as rQA part with a common ro 
for all systems and all incident energies. 

This formula, simplified for demonstration purposes, neglects other combinations 
leading to the production of K + as well as the use of Bose-Fermi distributions, which 
are all included in the computation. The corresponding formula for K~ production 
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is similar, but does not depend on (Xb- This point will become important later on. 

These formulae lead to a reduction of K + and K~ yields as compared to the 
numbers calculated without exact strangeness conservation [Q, [5j. Two extreme 
conditions can be seen from these equations. In the limit of a small number of 
strange particles, the additional term (due to the parameter V) leads to a linear rise of 
Mx+ /Ap ar t while M T / A part remains constant. This is in remarkable agreement with 
the experimental observations shown in Fig. [TJ. For very high temperatures or very 
large volumina, the terms in brackets approach unity (see Ref. [0) and the formulae 
coincide with the grand-canonical procedure. This is much better seen in the exact 
formulae using modified Bessel functions 0, 0, [ll[] . 

At low incident energies, the particle ratios (except tj/ttq) are well described using 
the canonical approach 14]. Surprisingly, even the measured K + / K~ ratio is described 
and this ratio does not depend on the choice of the volume term V. It should be noted 
that the statistical model uses normal masses of the particles while many transport 
calculation have to reduce the K~ mass (as expected for kaon in the nuclear 
medium) in order to describe the measured yields. It is therefore of interest to see 
how the results of the statistical model appear in a representation where the K + and 
K~ multiplicities are given as a function of ^/s — ^fsth as shown in Fig. ||. The choice 
of the x-axis follows the idea to correct for the different energy thresholds to produce 
K+ in NN collisions (y^ = 2.548 GeV) and K~ (</si/T = 2.87 GeV). It turned out 
that in this representation the measured yields of K + and K~ close to threshold in 
heavy ion collisions are about equal while they differ in pp collisions by factors of 10 

-loop. _ 

Figure |^ demonstrates that at values of yfs — y/sth less than zero, the excitation 
functions for K + and K~ cross leading to the observed equality of K + and K~ at 
SIS energies. The yields differ at AGS energies by a factor of five. The difference 
in the rise of the two excitation functions can be understood by the formulae given 
above. The one for K + production contains (E\ — /i^) while the other has E K + in 
the exponent of the second term. As these two values are different, the excitation 
functions, i.e. the variation with T, exhibit a different rise. 

Furthermore, the two formulae predict that the K + /K~ ratio for a given collision 
should not vary with centrality as V cancels in the ratio. Indeed, this has been observed 
in Au+Au/Pb+Pb collisions between 1.5 A-GeV and RHIC energies g| [jj, 0, ||, 
as shown in Fig. |^. This independence of centrality is most astonishing as one 
expects at low incident energies an influence of the different thresholds and the density 
variation with centrality. At 1.93 A-GeV the K + production is above and the K~ 
production below their respective NN thresholds. 
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Figure 2. Calculated K + /A part and 
K~ /Ap ar t ratios in the statistical model 
as a function of ^/s — y/sth for Ni+Ni 
collisions. The points are results for 
Ni+Ni collisions at SIS energies [P %. P 
and Au+Au at 10.2 A GeV (AGS) [fg\ 
At AGS energies the influence of the 
system mass is negligible. 
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Figure 3. The K+/K~ ratio appears 
to be constant of a function of centrality 
from SIS up to RHIC energies. The 
dotted lines represent the predictions of 
lie statistical model. Data from JlEj, 



Transport-model calculations show clearly that strangeness equilibration requires 
a time interval of 40 - 80 fm/c |2(], On the other hand, the statistical models 
assuming chemical equilibration are quite successful in describing the particle yields, 
including strange particles. 

In case of the K + production, no strong absorptive channel seems to be available 
which could lead to chemical equlibration. For K~ production the situation is quite 
different. At low incident energies strange quarks are found only in a few hadrons. The 
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Figure 4. Measured K /K + ratio, M(w° + n )/A par t and the double ratio 
([K-]/[K+])/([M(tt° + tt~)]/ A P art) as a function of A part ) both for Ni+Ni and 
Au+Au collisions at 1.5 A-GeV. Preliminary results!. 



s quark is essentially only in K + , while the s quark will be shared between K~ and 
A (or other hyperons). This sharing of the s quark might be in chemical equlibrium 
as the reactions 

ir° + A^p + K~ or 7r~ + A ^ n + K~ 

are strong and have only slightly negative Q-values of -176 MeV. 

The idea that the K~ yield is dominated by strangeness exchange via the tt~ + A 
channel has been suggested by p2[ and has been demonstrated quantitavely in a recent 
theoretical study |23|] . The direct K + K~ pair production via baryon-baryon collisions 
has negligible influence as these K~ are absorbed entirely. In these transport-model 
calculations the strangeness exchange is approaching equilibrium but does not fully 
reach it p3| . 

If these reactions are the dominating channel, the law of mass action might be 
applied giving for the respective concentrations pij 

M-[A] _ Jr 

7 ^ — /v. 

[K~] ■ N 

As the number of K~ relative to A is small, [A] can be approximated by [K + ] and 
rewriting gives 

|— | <XM{TT° +TT-)/ A par t. 

This relation also explains the measured constant ratio of K~ / K + with centrality 
(Fig. |§|) as the pion multiplicity does not vary with centrality. Figure |] demonstrates 
the constancy of the K~ / K + ratio and of the pion multiplicity with A par t for Ni+Ni 
and Au+Au collisions at 1.5 ^4-GeV Jl^, [l7|, It turns out that these ratios do 
not even depend on the choice of the collision system. The right part of this figure 
exhibits the double ratio ({K~]/[K + })/ ({M(tt° +ir~)]/A part ) which shows only a minor 
deviation from a horizontal line. This result can be taken as an argument that this 
specific channel might be not far from chemical equilibrium. 

Next wc test in Fig. || the validity of the law of mass action by plotting the 
K~/K + ratio as a function of the pion multiplicity M(tt° + n~)/A part at incident 
energies from SIS up to RHIC. At SIS and AGS energies the direct relation holds, 
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Figure 5. The K~/K~*~ ratio as a 
function of the pion multiplicity M(w~ + 
tt )/ Apart as a test of the law of mass 
0.1 1 10 action. Preliminary data. The dashed 

p. _ line shows the prediction of the statistical 

M(jI U +JI ) /A part model. 

i.e. the K~ / K + ratio rises with M(tt° + 7T~ ) /A part . At SPS and RHIC energies K 
are obviously produced by other channels, i.e. K + K~~ pair production. This change 
of the dominating channel is well reproduced by the statistical model (dashed line in 
Fig. |). 

2.2. Maximum relative strangeness content in heavy ion collisions around 30 A-GeV 

The experimental data from heavy ion collisions show that the K + /ir + ratio rises 
from SIS up to AGS but it is larger for AGS than at the highest CERN-SPS energies 

H n, ||, n EH and even at RHIC 11 as shown in Fi e- & 

This behavior is of particular interest as it could signal the appearance of 
new dynamics for strangeness production in high energy collisions. It was even 
conjectured Q that this property could indicate an energy threshold for quark-gluon 
plasma formation in relativistic heavy ion collisions. 

In the following we analyze the energy dependence of strange to non-strange 
particle ratios in the framework of a hadronic statistical model. In the whole 
energy range, the hadronic yields observed in heavy ion collisions resemble those of 
a population in chemical equilibrium along a unified freeze-out curve determined by 
the condition of fixed energy /particle ~ 1 GeV || providing a relation between the 
temperature T and the baryon chemical potential fis- As the beam energy increases 
T rises and /xs is slightly reduced. Above AGS energies, T exhibits only a moderate 
change and converges to its maximal value in the range of 160 to 180 MeV, while hb 
is strongly decreasing. 

Instead of studying the K + /ir + ratio we use the ratios of strange to non-strange 
particle multiplicities (Wroblewski factor) J3l| defined as 

2(ss) 
s (uu) + (do 1 ) 

where the quantities in angular brackets refer to the number of newly formed quark- 
antiquark pairs, i.e. it excludes all quarks that were present in the target and the 
projectile. 
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Figure 6. K + /ir + ratio obtained 
around midrapidity as a function of 
y/s from the various experiments. The 
dashed line shows the results of the sta- 
tistical model in complete equilibrium. 



Applying the statistical model to particle production in heavy ion collisions 
calls for the use of the canonical ensemble to treat the number of strange particles 
particularly for data in the energy range from SIS up to AGS j|, ^2) as mentioned 
before. The calculations for Au-Au and Pb-Pb collisions are performed using a 
canonical correlation volume defined above and given by a radius of ~ 7 fm determined 
in Q] . The quark content used in the Wroblewski factor is determined at the moment of 
chemical freeze-out, i.e. from the hadrons and especially, hadronic resonances, before 
they decay. This ratio is thus not an easily measurable observable unless one can 
reconstruct all resonances from the final-state particles. The results are shown in 
Fig. [?] as a function of yfs. 

The solid line (marked "sum" ) in Fig. ^describes the statistical- model calculations 
in complete equilibrium along the unified freeze-out curve |jj and with the energy- 
dependent parameters T and fiB- From Fig. |^ we conclude that around 30 A-GeV 
laboratory energy the relative strangeness content in heavy ion collisions reaches a 
clear and well pronounced maximum. The Wroblewski factor decreases towards higher 
incident energies and reaches a limiting value of about 0.43. For details see Ref. [^3| . 

The appearance of the maximum can be traced to the specific dependence of \ib 
and T on the beam energy. Figure^ shows values of constant X s in the T — hb plane. 
As expected X s rises with increasing T for fixed [Ib- Following the chemical freeze-out 
curve, shown as a full dashed line in Fig. g[ one can see that \ s rises quickly from 
SIS to AGS energies, then reaches a maximum around /j,b ~ 500 MeV and T w 130 
MeV. These freeze-out parameters correspond to 30 GeV laboratory energy. At higher 
incident energies the increase in T becomes negligible but \ib keeps on decreasing and 
as a consequence \ s also decreases. 

The importance of finite baryon density on the behavior of A s is demonstrated 
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Figure 7. Contributions to the Wrob- 
lewski factor A s (for definition see text) 
from strange baryons, strange mesons, 
and mesons with hidden strangeness. 
The sum of all contributions is given by 
the full line. 



Figure 8. Lines of constant Wroblewski 
factor A s (for definition see text) in the 
T — [ib plane (solid lines) together with 
the freeze-out curve (dashed line) 0. 



in Fig. ^ showing separately the contributions to (ss) coming from strange baryons, 
from strange mesons and from hidden strangeness, i.e. from hadrons like <f> and rj. 
As can be seen in Fig. the origin of the maximum in the Wroblewski ratio can 
be traced to the contribution of strange baryons. This channel dominates at low y^s 
and loses importance at high incident energies. Even strange mesons exhibit a broad 
maximum. This is due to the presence of associated production of e.g. kaons together 
with hyperons. 

The energy dependence of the K + /n + ratio measured at midrapidity is shown 
in Fig. |^. The model gives an excellent description of the data, showing a broad 
maximum at the same energy as the one seen in the Wroblewski factor. In general, of 
course, statistical-model calculations should be compared with 47r-integrated results 
since strangeness does not have to be conserved in a limited portion of phase space. 
A drop in this ratio for An yields has been reported from preliminary results of the 
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NA49 collaboration at 158 AGeV p8[ . This decrease is, however, not reproduced by 
the statistical model without further modifications, e.g. by introducing an additional 
parameter 7 S ~ 0.7 This point might be clearer when data at other beam energies 
will become available. 



3. Summary 

Strange particle production in heavy ion collisions close to threshold can be 
described by a statistical model in canonical formulation, i.e. using local strangeness 
conservation. This approach is able to explain several features of K + and K~ 
production at SIS energies. 

While for K + production it remains open whether and how chemical equlibrium 
can be reached, the situation for K~ is quite different. It is shown that the strangeness 
exchange process 7rA ^ N + K~ is the dominant channel for K~ production at SIS 
and likely also at AGS energies. This is demonstrated by applying the corresponding 
law of mass action. Theoretical studies confirm this interpretation. 

Using the energy dependence of the parameters T and /ib we have shown that 
the statistical-model description of relativistic heavy ion collisions predicts that the 
yields of strange to nonstrange particles reaches a well defined maximum near 30 GeV 
lab energy. It is demonstrated that this maximum is due to the specific shape of the 
freeze-out curve in the T — \ib plane. In particular, a very steep decrease of the baryon 
chemical potential with increasing energy causes a corresponding decline of relative 
strangeness content in systems created in heavy ion collisions above lab energies of 30 
GeV. The saturation in T, necessary for this result, might be connected to the fact 
that hadronic temperatures cannot exceed the critical temperature T c ~ 170 MeV for 
the phase transition to the QGP as found in solutions of QCD on the lattice. 
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